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The contribution from the R parity violating interaction, \[j^.LiQjD'f., in the associated produc- 
tion of a top quark (antiquark) with a charm antiquark (quark) is examined for high energy leptonic 
colliders. We concentrate on the reaction, l~ +1'^ -^ {tc) + (ct) -^ (blue) + (1)11/0), associated with the 
semileptonic top decay. A set of characteristic dynamical distributions for the signal events is eval- 
uated and the results contrasted against those from the standard model W-boson pair production 
background. The sensitivity to parameters (R parity violating coupling constants and down-squark 
mass) is studied at the energies of the CERN LEP-II collider and the future linear colliders. Next, 
p^ ■ we turn to a study of a CP-odd observable, associated with the top spin, which leads to an asym- 

f^ ' metry in the energy distribution of the emitted charged leptons for the pair of CP-conjugate final 

^\| , states, blfc and bWc. A non vanishing asymmetry arises from a CP-odd phase, embedded in the R 

parity violating coupling constants, through interference terms between the R parity violating am- 
plitudes at both the tree and loop levels. The one-loop amplitude is restricted to the contributions 
from vertex corrections to the photon and Z-boson exchange diagram. We predict unpolarized and 
polarized rate asymmetries of order O(10~^) — O(10~^). An order of magnitude enhancement may 
be possible, should the R parity violating coupling constants X[jf. exhibit a hierarchical structure in 
^vj . the quarks and leptons generation spaces. 
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^ . I. INTRODUCTION 

JTT^' The flavor non diagonal fcrmion-antifcrmion pair production, l^l^ — )■ fjfj', where J ^ J' are flavor labels, 

(-H I represents a class of reactions where the high energy colliders could contribute their own share in probing new physics 
incorporating flavor changing and/or CP violation effects. As is known, the standard model contributions here 
are known to be exceedingly small, whereas promising contributions are generally expected in the standard model 
^^ , extensions. (Consult ref. M for a survey of the literature.) Of special interest is the case where a top quark (antiquark) 
^ • is produced in association with a lighter (charm or up) antiquark (quark). The large top mass entails a top lifetime, 
5r ' Ttop — [1-56 GeV{ ^^^^^y )'^]^^, significantly shorter than the QCD hadronization time, 1/Aqcd, which simplifies 
the task of jet reconstruction. Q The top polarization effects also constitute a major attraction, ^-^j The large top 



mass entails a spin depolarization time of the top which is longer than its lifetime, Tdepoi = [l-7MeV{^^)] > Tt, 



topj 



thus providing an easy access to top polarization observables. Polarization studies for the top-antitop pair production 
reaction, in both production and decay, have been actively pursued in recent years. [pl^lOl (An extensive literature 
can be consulted from these references.) 

It appears worthwhile applying similar ideas to the flavor non diagonal fermion pair production process involving 
a single top production. This reaction has motivated several theoretical studies aimed at both leptonic (l~l^, ej 
and 77) and hadronic (pp, pp) colliders. Exploratory theoretical studies have been pursued at an implicit level, via 
the consideration of higher dimension contact interactions [0-n3[ , and at an explicit level, via the consideration of 
mechanisms involving leptoquarks, p4[ an extended Higgs doublet sector, |l5^1q supersymmetry based on the minimal 
supersymmetric standard model with an approximately broken R parity, [|1[|17|-20|] quark flavor mixing, pl| standard 
model loops and four matter generations, |23,E3] or higher order standard model processes with multiparticle final 
states, l^l^ — > tcvP. pi] A survey of the current studies is provided in ref. [p!2[. 
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In this work, pursuing an effort started in our previous paper, ^ we consider a test of the R parity violating 
(RPV) interactions aimed at the top-charm associated production. Our study will focus on the contributions to the 
process, l^l^ -^ (tc) + (ic), arising at the tree level from the trilinear RPV interactions, \[jj^LiQjD%, via a dkR squark 
exchange. We examine the signal associated with the (electron and muon) charged semileptonic decay channel of 
the top, t -^ bW^ -^ bl^v. The final states, {bl^vc) + (bl~vc), [I = e, fx] consist of an isolated energetic charged 
lepton, accompanied by a pair of b and c quark hadronic jets and missing energy. The standard model background 
may arise from the W-boson pair production reaction, l'^l~ — > W~^W~ , and possibly, in the case of an imperfect b 
quark tagging, from the b — b quark pair production reaction, I'^l^ —> bb, followed by a semileptonic decay of one of 
the b quarks, b -^ cl^D. 

The present work consists of two main parts. In the first part, we discuss the signal associated with the top 
semileptonic decay channel. We evaluate a set of characteristic dynamical distributions for the signal and for the 
standard model background and obtain predictions for the effective rates based on a judicious choice of selection cuts 
on the final state kinematical variables. Our discussion will develop along similar lines as in a recent work of Han 
and Hewett, |12| which was focused on the contributions initiated by the dimension, I? = 6, four fields couplings of 
the Z-boson with fermion pairs and the neutral Higgs boson. In the second part of the paper, we examine a specific 
CP-odd top polarization observable which corresponds to an asymmetry in the energy distribution of the final state 
charged lepton with respect to the sign of its electric charge. 

The contents of the paper are organised into 3 sections. In Section^, we focus on the total and partial semileptonic 
decay rates for both the signal and standard model background, allowing for the case of an imperfect b quark tagging. 
We discuss the constraints from the indirect bounds on the RPV coupling constants, study the dependence of rates 
on the down-squark mass parameter and evaluate a set relevant dynamical distributions that are of use in devising 
an appropriate set of selection cuts. In Section HI, we discuss a test of CP violation involving top polarization 
effects. The CP violating observable arises through interference terms between the tree and one-loop contributions 
to the amplitude and a CP-odd phase which is embodied in the RPV coupling constants. Following an approach 
similar to one used in earlier proposals, |gj7[] we describe the top production and decay by means of a factorization 
approximation and examine the induced charge asymmetry in the energy distribution of the final state charged leptons. 
The production amplitudes are evaluated in the helicity basis. Our main conclusions are summarized in Section IV. 



II. TOP-CHARM ASSOCIATED PRODUCTION 



A. Integrated rates 



In a l^l'^ collision, the tree level transition amplitude for single top production, as initiated by the RPV interactions, 
X'^Aj.LiQjD'f., proceeds via the u-channel exchange of a right-handed down-squark, dkn, as represented in Fig|l|. By use 
of a Fierz ordering identity, the transition amplitude for the flavor non diagonal production of an up quark-antiquark 
pair, l^ (k) + l'^{k') -^ uj{p) + uji{p'), can be written in the form of a Lorentz covariant vectorial coupling. 
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FIG. 1. Feynman diagram for the tree level amplitude of the process, l^l^ -^ ct -^ cbl^v. 



We shall specialize henceforth to the case of electron-positron colliders, corresponding to the choice I = 1 for the 
generation index. The squared amplitude, summed over the initial and final fermion spins, reads: |lj] 
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The production rate for unpolarized initial Icptons, integrated over the scattering angle in the interval, < | cos 0\ < Xc, 
is given by the analytic formula, 

- - ^-\^'iJ'k^iJk\\^_ _ „^) + (2^2 _ ^2 _ ^2,) In |: 
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647rs' 
- {rh' - m]){rh' - m^K ^ ^)], (3) 

where, u± = ?7ij — ^/s{Ep -izpxc)- For the top-charm associated production case, in the limit, mj = rrit >> mji — rric, 
one has, u+ ~ m^ — s, u_ ~ 0. For fully polarized initial beams, since the RPV amplitude selects a single helicity 
configuration for the initial state leptons, I^^ri (l^ft handed l~ and right handed l^) the corresponding polarized rate 
would be still described by the same formula as above, only with an extra enhancement factor of 4. The predicted rates 
for tc production are controlled by quadratic products of the RPV couphng constants, ^i3kK2k^ [^ = 1' 2, 3] and the 
squark mass, toj , denoted for short as, tti. Allowing for the existence in the RPV interactions of an up-quark flavor 
mixing, such as would be induced by the transformation from flavor to mass basis, one may express the amplitude in 
terms of a single RPV coupling constant and the CKM matrix, V, by rewriting the coupling constant dependence as, 
'^I2fe'^''^i3fe ~^ '^iMfc'^iM'fc(^^)M'2(^^)M3> ^^^ Selecting the maximal contribution associated with the configurations, 
M = M' = 2 or 3. This yields the order of magnitude estimate, X'i%kK2k ^ l^i2feP(^''')22(V^^)23 ~ '^\K2k\'^^'^ or 
2|A'j^3j,pA^, respectively, where, A « sin^c ~ 0.22, denotes the Cabibbo angle parameter. 

We pause briefly to recall the current bounds on the RPV coupling constants of interest in the present study. 
p5| The relevant single coupling constant bounds are, A'j^j/c < 4. x 10^^, A'^^g^ < 0.37 (charged current universality); 
X'lji < 3 X 10~^ (atomic physics parity violation); A'j2fe < 0-3 ~ 0.4, A'^gj, < 0.3 — 0.6 (neutral current universality); 
and Ai22 < 7. x 10^^, A'^gg < 3.5 x 10^'' (neutrino Majorana mass). |^ The superpartners scalar particles masses are 
set at 100 GeV . Unless otherwise stated, all the dummy flavor indices for quarks and leptons are understood to run 
over the three generations. Using the above results for individual coupling constants bounds, we may deduce for the 
following upper bounds on the relevant quadratic products, ^ KskKik < [0(10"^), 0(10"^), 0(10"'')], [k = 1, 2, 3]. 
The indirect quadratic products bounds, X^j,^X[,^,. < 1.1 x 10~^, X[j,^X[,^^ < 1.1 x 10~^, [i' = 1,2] (J5 -^ Xqvv) are 
roughly comparable to these single coupling constants bounds. We also note that using the CKM flavor mixing along 
with a single dominant coupling constant in the current basis, as described at the end of the previous paragraph, 
may not be especially beneficial in avoiding the above stronger pair product bound. The bound on the corresponding 
coupling constant factor, 2]A'j^3j,]^A^ < 0(10^^), is competitive for the generation indices, fc = 1,3. 

Numerical results for the integrated rates have already been reported in previous works ]^y,y . Setting the relevant 
RPV coupling constants at the reference value. A' = 0.1, one predicts rates of order 1 — 10 /6, for m = O(IOO) GeV. 
As the center of mass energy varies in the interval, -y/s = 192 — 1000 GeV, the rates rise sharply from threshold, 
reaching smoothly a plateau around ^/s ~ 400 GeV . This contrasts with the predictions from gauge boson mediated 
higher dimension interactions ]|l3] where the rise of the rates with incident energy is a more gradual one. The rates 
are also found to have a strong dependence on m, which weakens for increasing center of mass energies. One may 
roughly parametrize the dependence on s and tti by the approximate scaling law, a « {jrkiY { ^^^ ^'^^ Y^''\ where the 
power exponent is a fastly decreasing function of energy, taking the approximate values, x{s) w [3.65, 1.86, 0.94], at, 
Vi= [0.192, 0.5, 1.] TeV. 

Although the predicted rates seem to be severely constrained by the above indirect bounds, one could envisage an 
optimistic scenario where the supersymmetry decoupling limit, fh —f oo, is realized with fixed values for the products, 
A^ j,(i22^^£K) ~ 0.1, consistently with the current indirect bounds. The results obtained with this prescription are 

displayed in Figure ||. The integrated rates now depend on m as, u ex { ^^^ ^''^ )~^~^^'--^\ which leads at high energies 
to an enhancement by up to three orders of magnitudes, compared to the case where the RPV coupling constants 
are taken independent of rh. The initial energy of LEP-II falls right in the regime where the cross section is sharply 
rising with increasing initial energy. The decrease with increasing m is stronger at LEP-II energies than at the future 
linear colliders energies. Note that at the largest values of the superpartner mass, ttt, ~ 1. TeV, the RPV coupling 
constants in our prescription enter a strong coupling regime (A' = 0(1)) and it is not clear then whether the tree level 
prediction makes sense. 
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FIG. 2. The total integrated rate for the RPV induced reaction, I'^l —>■ (tc) + (tc), setting the values of the relevant coupling 
constants as, \[2k = -^'isfc ~ 0.1 (m/100 GeV), is plotted in window (A) as a function of center of mass energy, s^^^ for fixed 
down squark mass, m = [100, 200, 500, 1000] GeV and in window {B) as a function of m for fixed s5 = [192, 500, 1000] GeV. 
We integrate over an interval of the scattering angle, < | cos6\ < 0.9848, corresponding to an opening angle with respect to 
the beams axis larger than 10°. 

Next, we consider the process incorporating the top semileptonic decay, as pictured by the Feynman diagram 
shown in FigJ^. We assume that the top decay is dominated by the electroweak semileptonic decay channel, with 
branching fraction, B{t ^ b + W~^) « 1. We also include the pair of CP-conjugate final states, tc and ct production, 
which multiphes the rate by a factor of 2. Note, however, that we restrict ourselves to the uj' = c charm quark 
mode only. The numerical results for rates, including a branching fraction factor of | (experimental value, 21.1%) 
to account for the W -^ Iv, [I — e, /i] decay channels, are displayed in Table II A. We also show the standard 
model background rate from the W^-boson pair production, ?^/+ -^ [W^W^) -^ ( l^vuidj) + {l~DdjUi), with one 
VF-boson decay leptonically and the other hadronically, where z, j are generation indices. The irreducible background 
from, W^ ~* cb or W~^ — *■ &c, is strongly suppressed, due to the small branching factor, given approximately by, 
0.32|Vcb|^ ~ 5 10^''. It is safer, however, to allow for the possibility where the light quark hadronic jets could be 
misidentified as b quark hadronic jets. Accounting for the leptonic decay for one of the M^-boson and the hadronic 
decay for the other VF-boson, introduces for the total rate, which includes all the subprocesses, the b ranching fraction 
factor, 2 x (21.1 ± 0.64)% x (67.8 ± 1.0)% = 0.286 ± 0.024. Our numerical results in Table |l^ for the standard 
model background rates are in qualitative agreement with those quoted {a = [2252, 864] fb at s^ — [0.5, 1.] TeV) by 
Han and Hewett. [|2| One should be aware of the existence of large loop corrections to the W~^W~ production rate, 
especially at high energies. The predictions including the electroweak and QCD standard model one-loop contributions 
read, (28| cr = [4624, 1647, 596] fb at s^ = [0.192, 0.5, 1.] TeV. We conclude therefore that our use of the tree level 
predictions for the W^W^ background overestimates the true cross sections by [9%, 20%, 32%] at the three indicated 
energies. 



TABLE I. Production rates for the top-charm production signal and the W-boson pair production background. The line 
entries give successively the total integrated rate for the reaction, l^l" -^ (tc) + (cf), using, A' — 0.1, m = 100 GeV, the rate 
for signal events, (bli/c) + (bWc), associated with the top semileptonic decay, the W-boson pair production background rate, 
l~l'^ —f W^W' —f {I'^vUidj) + {l''vdjUi), and the corresponding cut signal and background rates, as obtained by applying the 
selection cuts quoted in the text. The results include the first two generations of charged leptons, I — e, /i. 



Energy(TeV) 


0.192 


0.5 


1.0 


Total rate a{fb) 


4.099 


4.291 


1.148 


Signal (fb) 


0.68 


0.91 


0.24 


W+W- Background (fb) 


5076 


2080 


876 


Cut Signal (fb) 


0.54 


0.74 


0.21 


W+W~ Cut Background (fb) 


17. 


5.0 


2.6 



Let us discuss briefly other possible sources of background. The next important contribution is that arising from the 
non resonant W-boson virtual propagation in the amplitude with the intermediate W~^W~ bosons branching into four 
fermions [Ivqq'). This could be possibly estimated by subtracting the resonant contribution from the total background 
cross sections, weighted by the suitable branching factors, as independently evaluated by numerical methods in the 
literature. The results for the integrated total cross sections, l~l'^ -^ (4/) + (4/ + 7), |29] including the initial state 
radiation and Coulomb corrections, indicate that the off-shell contributions amount to a small relative correction lower 
than 0(10%). Alternatively, one may consider, after reconstructing the neutrino momentum from the missing energy, 
a procedure to impose suitable cuts on the hW invariant mass, aimed at suppressing the non resonant production 
background. 

One other potentially important background is that arising from the h — b quark pair production reaction, Z+Z~ — > 
7*/Z -^ bb ^ b{cl~i>) + b{cl^v). iQ The numerically derived predictions for the rates, as obtained by means of 
the PYTHIA generator, are: a = [1.631 10^, 2.12 10^, 5.35 10^] /6, at s* = [0.192, 0.5, l.]TeV. It would appear 
desirable, in view of these large predicted rates, to eliminate this background by performing a double b quark tagging 
analysis on the events sample. This can be performed at a reasonably low cost, given that the detection efficiency of 
b quark jets is currently set at 50%. If one performs a single b quark tagging, the rates for the corresponding events, 
l^l~ -^ 1* /Z -^ bb -^ {cl~D)b, are reduced by a branching fraction, B{b -^ civ) « 10%, but this is compensated 
by the probability of misidentifying a light quark jet as a 6 quark jet, which lies at the small value of 0.4% with the 
current silicon vertex techniques. If no b quark tagging is performed at all, then the above large rates may make it 
necessary to resort to an analysis of isolation cuts of the type to be discussed in the next subsection. 

B. Distributions for the semileptonic top decay events 

In order to separate the signal from background, we consider the same set of characteristic final state kinematical 
variables as proposed in the study by Han and Hewett. []13| These are the maximum and minimum energy of the two 
jets, E"^^, EY^" , the dijet invariant mass, Mjj, the charged lepton energy, Ei, and rapidity, yi ~ j log g'_^'" ■ The 

distributions in these 6 variables for the signal and background, at a center of mass energy, y^ — 0.5TeT^, are plotted 
in Fig.|[ These numerical results were obtained by means of the PYTHIA [Q event generator. One notices marked 
differences between signal and background. The maximum jet energy distribution is uniformly distributed for the 
background but sharply peaked for the signal, where the peak position is determined by the top mass and the incident 
energy as, m^ = (m^ — s + 2y^Ep). The minimum jet energy is uniformly distributed for both signal and background, 
but happily the corresponding intervals are very partially overlapping. The signal events rapidity distributions for 
the maximum energy jet are more central for signal than background. A similar trend holds for the lepton rapidity 
distributions. The dijet invariant mass is a most significant variable in discriminating against the background due 
to its pronounced peak at the W mass. For the signal, the dijet invariant mass is uniformly spread out. Although 
we do not show here the distributions for the top mass reconstruction, this also features a strong contrast between 
a strongly peaked signal and a uniform background. The lepton energy distributions for the signal and background 
are peaked at the opposite low and high energy ends of the physical interval, respectively. This is a familiar effect 
associated with the correlation between the W-boson spin polarization, which is predominantly longitudinal in the 
top decay and transverse in the direct W-boson decay, and the velocity of the emitted charged lepton. In the signal 
decay amplitude, t —^ blv, the fact that the left handed b-quark must carry the top polarization, forces the lepton 
to travel with opposite velocity to that of top. In the background decay amplitude, W~ -^ W, the charged lepton is 
emitted with a velocity pointing in the same direction as that of the W-boson. Thus, the Lorentz boost effects on the 
emitted charged leptons act in opposite ways for the signal and background events. 

While the above distinctive features between signal and background events get further pronounced with increasing 
center of mass energy, opposite trends occur as the initial energy is lowered. The distributions at the LEP-II center 
of mass energy, ^/s = 0.192 TeV, are plotted in FigJj. At this energy, the monovalued distribution for the signal jet, 
which is now the softer lower energy jet, is still well separated from the corresponding background jet distribution. So, 
this variable, along with the dijet invariant mass stand up as useful discrimation tests for the signal. By contrast, the 
energy and rapidity distributions for the maximum signal jet may not be easily distinguished from the background. 
Similarly, the lepton energy distributions in the signal and background are overlapping due to the small Lorentz boost 
effect. 

The distributions obtained with the RPV interactions are rather similar to those found with the higher dimension 
operator mechanism. [jl3| This is due to the formal structure of the RPV amplitude, involving an effective u-channel 
vector particle exchange. In fact, the selection cuts proposed by Han and Hewett |12| appear to be quite appropriate 
also in the RPV case, and, for convenience, we recapitulate below the cut conditions used to characterize the selected 
events. 



£;j°'" < 20, £;f ^'' > 60, El > 0, Sjj > 10, St < 5, [Vi = 192] 

E'-°"' > 20, £;f ^'' > 200, El < 150, 6jj > 10, (5f < 40, [^s = 500] 

£;j°™ > 20, Ej'^'' > 460, El < 350, <5jj > 10, St < 100. [^/i = 1000] 

The above listed variables correspond to the minimum and maximum energy of the two jets, £''°"', E"^'', the 
charged lepton energy, Ei, the distance between the dijet invariant mass and W-boson mass, djj — \Mjj — rnwl, 
the distance of the reconstructed top mass to the true mass, St = \rn'['^^°'^^*' — rat\. The assigned numerical values 
are all expressed in GeV units. Besides the above cuts, we also impose the usual detection cuts on energies and 
rapidities, Ej,i > 10 GeV, {ijj^il < 2, aimed at removing the particles travelling too close to the beam pipe. We 
allow for the detection efficiency of the particle energies only in an approximate way, namely, by accounting for the 
following approximate uncertainties, AE/E = 40%, 10%, on the jets and lepton energies, respectively, at the level 
of imposing the above selection cuts, rather than by the usual procedure of performing a Gaussian smearing of the 
particle energies. 

The numerically evaluated efficiencies on the signal and background events are, eg ~ 0.8, es ~ 3.10^'^, with a 
very weak dependence on the center of mass energy and, for the signal, a weak dependence on the mass parameter 
TO, which was set at m = 100 GeV in the numerical simulations. After applying the cuts, the background rates are, 
ases ^ [17., 5., 2.] fb, and the signal rates, as^s = [0.68, 0.74, 0.21] fb, for ^/i = [192., 500., 1000] GeV. The results 



for the cut signal and background rates, as given in Table [1 A , show that the background is very significantly reduced 



by the cuts. The situation is clearly far more favorable for the future linear colliders than for LEP-II. Nevertheless, 
the number of surviving signal events is still one order of magnitude below that of the surviving background, so that 
the option of cutting down the background by means of a 6 quark tagging procedure is to be preferred since the 
ensuing reduction would be much more drastic. An integrated luminosity oi C = 100 fb^^ would lead to a number of 
signal events, iK2kK3k/'^'^~'^) ^ 0(30). 

We have also performed an indicative event generator study of the background, l^l^ -^ bb ^ I + hadrons, 
restricting consideration to the emitted charged leptons only. A jet reconstruction of the partonic level distributions 
is a task beyond the scope of the present work. We focus on the first charged lepton emitted during the semileptonic 
decays of the produced B, B mesons, since this carries the largest velocity. As seen on Fig.pl the distribution for 
the first emitted charged lepton energy is peaked at low energies. One expects that the most energetic lepton is that 
produced in the semileptonic decays of the B mesons. The rapidity distribution is less central than for the signal 
and nearly overlaps with that of the W^W~ background. Therefore, imposing the additional lower bound cut on 
the lepton energy, say, at Ei > 20 GeV, for 52 = 500GeV, should be sufficient to appreciably suppress the b — b 
background without much affecting the signal. 

We may infer the reach with respect to the free parameters by evaluating the statistical significance ratio for a 
discovery, as defined by, a = -y===, S — asC, B = (tbC, where L denotes the integrated luminosity. Setting this at 

the value, a — Z, corresponding to a 95% confidence level, one deduces a dependence of the RPV coupling constant 
as a function of the superpartner mass parameters for a fixed initial energy and integrated luminosity. The sensitivity 
reach contour plot for the relevant parameters. A' A' = ^'i2k^vik ^^^ "^ ~ "^d ' ^^ shown in FigS We note that the 
sensitivity limit on the product of coupling constants, A' A', scales with the luminosity approximately as, 1/\AC. While 
the reach on the RPV coupling constants products, A'j^j/c-^ia/c ^ O(10~^), lies well above the current indirect bounds, 
this covers a wide interval of the down squark mass which extends out to ITeV. To compare with analogous collider 
physics processes, we note that while the flavor diagonal fermion pair production reactions, e~e"*" — > fjfj, may have 
a higher sensitivity reach, these are limited to information on the single coupling constants, \[jf. pM. The special 
reaction, e^e+ — > bb, proceeding via a sneutrino s-channel resonance, may probe quadratic products such as, AisiAggg^ 
| p2| or AiaiAj]^]^ [pSJ at levels of O(10~'^), but this is subject to the existence of a wide sneutrino resonance. The tb 
associated production at the hadronic Fermilab Tevatron [0,0 and the Cern LHC [Q colliders can be initiated 
via a charged slepton ej^ s-channel exchange. The sensitivity reach on the linear combination of quadratic coupling 
constants products, X'iiiX'iss, is of order 10~^ — 10^^. This information should prove complementary to that supphed 
by our study aimed at the leptonic colliders. To conclude this brief comparison, we observe that the information 
provided by the single top production reaction appears to be rather unique in view of the very characteristic signature 
of the associated events. 
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FIG. 3. Normalized dynamical distributions associated with the signal events, l^l -^ (tc) + (tc) -^ (W^vc) + {hi i^c), 
(dashed line) at m = 100 GeV, and the background events, I'^l^ -^ W^W~ —> {l^vqq) + [l^vqq'), (continuous line) at a 
center of mass energy, s^ = 500 GeV . The kinematical variables in the histograms, from left to right and up to down, are the 
jets maximum energy, the jets minimum energy, the rapidity for the highest energy jet, the di-jet invariant mass, the charged 
lepton energy and the charged lepton rapidity. The charged lepton energy and rapidity distributions are also plotted for the 
h ~b background production events, I'^l^ —* bb -^ l^ + hadrons (dotted line). 
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top — charm versus W* W at 1 92 GeV 
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FIG. 4. Same distributions as in Fig.H, at a center of mass energy s^ — 192 GeV. 
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FIG. 5. Sensitivity reach plot for the RPV coupling constants product, Ai2feAi3s./0.01, as a function of the down squark 
mass, m, for fixed center of mass energy, s' — [192. ,500, 1000] GeV, and corresponding fixed integrated luminosity, 
L — [2., 100., 100.] fb^^, using an acceptancy for the background, es = 3 10^^, and an acceptancy for the signal, es = 0.8, 
assumed to be independent of rh. 

III. TOP POLARIZATION OBSERVABLES AND A TEST OF CP VIOLATION 

Should single top production become experimentally observable in the future, an important next step to take is in 
examining top polarization observables. In this section, we present an approximate study for the top semileptonic 
decay signal in top-charm associated production aiming at a test of CP violation. We exploit an idea which was 
developed in early studies oi t — i production. |g,^ Interesting extensions are currently pursued ||8|-pi|. The basic 
observation is that any CP-odd quantity depending on the top polarization, such as the difference of rates between 
the pair of CP conjugate reactions, a{l^l^ -^ t^c) — a{l^l^ —> inc), can become observable by analyzing the 
top polarization through the kinematical distributions of its emitted decay (6 quark or charged lepton) products. An 
especially interesting observable is the charged lepton energy distribution for a polarized top. Any finite contribution to 
the CP-odd observables must arise through an interference term involving imaginary parts of loop and tree amplitudes 
factors, the loop amplitude factor bringing a CP-even final state interaction complex phase with the CP-odd relative 
complex phase arising from the coupling constants in the product of loop and tree amplitudes. 



A. Helicity basis amplitudes 

Building on our previous work, lit] we shall combine the tree level RPV induced amplitude discussed in Section 
II with the one-loop RPV induced amplitude associated to the photon and Z-boson exchange diagrams, restricting 
ourselves to the vertex corrections in the electroweak neutral current vertices, jf.j{p)fj'{p') and Zfj{p)fj'{p'). The 
Z-boson vertex admits the general Lorentz covariant decomposition. 
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-.JJ' 



2 cos 9w 



{Z). 



^.IJ' I 



U.I' 



J. J' I 



{Z)=^Mr {f)PL+A''^' {I)Pr) + 



1 



-<y^.u{p + Pr{^a"' +I5d''''), 



mj + mj 
where the vectorial vertex functions, ^^''r = A'l'^j^\tree + ^l'^rIioop, have a tree level contribution given by, A'l'^j^ 



(4) 



L,R.\tree 



nL,R 



5jj'aL,R{f), aL,Rif) = 2T3 ' — 2Q(/) sin 6w, and the tensorial vertex functions, a , d , are associated with the 
anomalous transition magnetic moment and the CP-odd, P-odd electric transition dipole moment, respectively. An 



analogous decomposition applies for the photon, J^ — _ g'^"^^w pJJ ^-^^^ with aL,Rif) = 2Q{f), determined by the 
electric charge Q{f). It is convenient to work with the Z fjfj' vertex in the alternate Lorentz covariant decomposition, 
r^'' {Z) = jf_i{A — B'j^) + |(p — p')fj.{C — 2575), where the vertex functions. A, B, C, V, (omitting the up quarks 
generation indices J, J' for convenience) are related to the previously defined vectorial and tensorial ones, eq. (]^, as, 

A ^ liAf'if) + A'/ if)) + a-''^\ B ^ \{Ai'\f) - Ak^'if)) + ''-ll-I!^^d"\ 

C = ^ a''', V = ? id'"'. (5) 

mj + mji mj + mj' 

The one-loop Z-boson exchange amplitude may then be written in the form, 

2 



Mi'''{Z)= (2^^) v{k',fi')j,(^a{eL)PL + a{en)PRy{k,fi) 



s 



^1 



X u{p,X)[r{A-Bj,) + ^ip-pr{C-V^5)]v{0,y)- (6) 

Combining the above loop amplitude with the RPV tree amplitude, eq. (O]), which we rewrite as, 



(2 \ ' 
u-mj_) 



M/'^ =7^^;7^(l-75)7.jI7''(l-75>, 7^^- ^,y^ ^f \ , (7) 



'J-kR 



one obtains, 

MJJ' ^ M/^' + Mf-''{Z) = [{Ga+A + 7^)(7^)(7^) - {Ga+B + n){j^){ri5) 
- {Ga-A + 7^)(7^75)(7") + {Ga-B + 7^)(7;.75)(7^75) 
+ ^(P -p')^[Ga+C(7^)(l) - Ga+P(7^)(75) - Ga-C(7;.75)(l) + Ga-V{^^^,)ij,)], (8) 

where, a — ^(aL(e)±afl(e)), and we have omitted writing the contractions of the Dirac spinors indices for the initial 
and final fermions, respectively. The photon exchange contribution can be incorporated by treating the parameters 
a as operators acting on the vertex functions. A, B, C, D, by means of the formal substitutions, 

J J' I 



Xe)±aR{e) / Aj^' jf) ± A^/ jf) ^ ^jA ^ ^^ /2Q{f)\ / AJ;'' jf) ± A],'' jf) 



Gz = (^^V)' 4-^' ^-^ = (^^)^i. (9) 

2 COS uw s — m^ + imzi^ z Is 

Analogous formulas to the above ones hold for the other products, Ga^B, Ga^C, Ga^V. We have labelled the vertex 
functions for the photon current by the suffix 7. The formulas expressing the RPV one- loop contributions to the vertex 
functions arc provided in Appendix |A|, quoting the results derived in our previous work. ||l|| The amplitude M'''^ in 
eq.(]q) may be viewed as a 4 x 4 matrix in the fermions polarization space, (/(p, A)/j'(p', A')|Af |Z+(fc', /i')/~(fc,/i)). 
The various products in eq. (g|) for the matrix elements with respect to the two pairs of Dirac spinors separate into 
8 distinct terms. The calculation of the helicity amplitudes is most conveniently performed with the help of the 
'Mathematica' package. Of the 16 configurations only the 8 helicity off diagonal configurations in the initial fermions 
are non vanishing. The explicit formulas for the helicity amplitudes are provided in Appendix El 

B. Charged lepton energy distribution 

The differential cross section for top production and decay is described in the factorization approximation. Ignoring 
the spin correlations, which corresponds to dropping the spin non diagonal contributions between the production and 
decay stages, yields: 

■'''-i2lWv/-''°°''')i:i"^.-^r'^-'>-')iV''''' b^-^;l.,,.,r.p -''-'' 

"^r* = o «» IZ„|Jfe.('A' ^ bl*v)\'dE;dEl (10) 

(iTTJ-^OTTlt -^ — 'A' 

10 



The production amplitude is denoted, Mprod, the top decay amphtude, Mdec, and A, A' = ±1, are polarization labels, 
which will also be written for short as, ± . We shall assume a narrow resonance approximation for the top propagator, 
|p^ — m^+ imtTt\~'^ -^ ^\ 5{p^ — m^ ). For the energies of interest, all the leptons and quarks, with the exception of 
the top, may be treated as massless. Two frames of interest are the laboratory {l~l~^) rest frame and the top rest frame. 
The letters denoting momentun variables in the l~l'^ center of mass (laboratonO frame are distinguished from those 
in the top rest frame by the addition of a star. Standard kinematical methods |3J| can be used to transform variables 
between these frames. Exploiting the rotational invariance, one may conveniently choose to work in the spatial frame 
where the top momentum lies in the xOz plane (0, = 0) and the charged lepton points in an arbitrary direction 
described by the spherical angles, Oi,(j>i. The relations between angles may be obtained by use of the spherical triangle 
identities, for example, the angle between lepton and top reads, cos6';j = cos 6*; 0036* + sm9i sin9 cos (j>i. The Lorentz 
boost from the top rest frame to the laboratory frame, involves a velocity parameter, v = p/Ep, j3 = p/Ep, 7 = 
(1 — /3^)~^" = Ep/rrit, and yields for the charged lepton momentum four vector and polar angle relative to the top 

momentum, E^ = jiEi - v ■ h), k^^ =ki +^v{^ - Ei), cos 0* = j^^, ■ 

The top differential semileptonic decay rate has been thoroughly studied in the literature. p5| One representation 
convenient for our purposes is the double differential rate with respect to the final charged lepton energy, E^ , and 
the final lepton and neutrino invariant mass squared, W^ = {ki + kj/)^- The result for the unpolarized rate carries no 
dependence on the scattering angles and reads, quoting from ref. [ p5[ , 

,^ NiGjpm^ f xi{xM-xi) 
"^t = — --r^ — dxi / dy 



^ NiG%ml 2 xijXM - xi) ^^^^_^ J^xiJxm - xi) ^^^ 

167r3 mt 7^ {l + j^){l-xi)-(,xi{xM-xi) 

The kinematical variables for the emitted charged lepton and neutrino are defined as, xi — 2Ei/mt, y = 
W'^ /ml, \W = ki + k^] with the bounds, < x; < xm, < y < ^ i^f^^"^ ^^^ ^^ employ the following nota- 
tions, Ni for the number of light lepton flavors, 7 = Tw/mw, S. = rnl/rriy^,, xm = 1 — e^, e = rnb/nit, tan"-'^ A = 
Artan\A\ +'k9{—A). Recall that the number of light lepton flavors, Ni, is set to iV; = 2 in our analysis. A useful trick 
to obtain the distribution with respect to the laboratory frame lepton energy, Ei, is to choose the top momentum 
along the Oz axis fixed frame and introduce the top rest frame electron energy by means of the change of variable, 
{Ei,cos9f) — > {Ei,Ef), associated with the Lorentz boost between the top rest frame and the laboratory frame. 
El = jEI{l + l3cos9l). The result reads, 

dTt f+\ d^Vt^ 2 /■< dxi d^Tt ,^^, 

= / d cos 9i — / , (12) 

dEi J_i dEidcos9i mtjP J^- xi dxidcos9i 

where the integration interval over xi is bounded at, x^ — — — nVgy. 

C. Top polarization observables 

An essential use will be made of the factorization property of the double differential distribution for the top decay 
semileptonic rate with respect to the emitted lepton energy and angle relative to the top spin polarization vector. 
This distribution is described at the tree level as, ^lE'^dcosib ~ jgr -'-+^°'^'^' ^ where, cosi/"; = ^s(p) • fc;, is the angle 
between the lepton momentum and the top spin polarization vector, s^{p), in the top rest frame. Equivalently, 
dx d cos e* ~ ^ ^'^'^T dmst* ' Ull ^® ^^ tums out, this representation remains vahd to a good approximation when 
one-loop QCD corrections are included. [Q We choose to describe the top polarization in the spin helicity formalism, 
using techniques familiar from previous works. |5|,|37| The definition for the helicity basis Dirac spinors is provided in 
Appendix H. Since the polarization axis coincides then with the top momentum, the dependence on ipi can also be 
simply rewritten as, (l-|-cos'0/)/2 = (l-|-Acos0*)/2, such that, A = [—1, -1-1], correspond to [L,R] hehcity, respectively. 

The helicity amplitudes associated to the pair of CP-conjugate processes are related by the action of CP as, 
< /a/a'I-^^IC'^^ >^'< f'_\'f-\\M\l^^lZij^i > ■ Unlike the process, l'^l~ -^ ti, where both the initial and final states 
are self-conjugate under CP, here only the initial state is self-conjugate, while the action of CP relates the different 
final states, tc and ci. Let us express the amplitudes for the pair of CP-conjugate processes as sums of tree and 
loop terms, M"^'^ = a^ -\- ^^ bafa{s + ie), M"^'^ = a^ + ^^ b*^fa{s + ie), where the loop terms, bafa{s + ie), are 
linear combinations with real coefficients of the vertex functions, Aj^'^ , Aj^ , a^^•^ , id"^'^ , with the energy dependent 
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complex functions, /a(s + ie), representing the factors in loop amplitudes which include the absorptive parts. In terms 
of these notations, a CP asymmetry associated with the difference of rates for the pair of CP-conjugate processes in 
some given CP-conjugate configurations of the particles polarizations, can be written schematically as, 

I < AA'|M|/iV > P - I < -A' - A|Af I - ^ - ^' > |2 ex ^ /m(ao6*)/m(/„(s + ie)) 

a 

- Y, lMbc.b*^')lMfa{s + te)f*,{s + ie)). (13) 

a<a' 

Thus, the necessary conditions for a non vanishing polarized asymmetry to arise from the tree-loop interference term 
are a relative complex CP-odd phase between the tree and loop coupling constants and an absorptive part from the 
loop terms. The angle integrated production rates for the CP-conjugate reactions, l^l~ -^ tc and l~^l~ -^ ci, for the 
case of polarized top and antitop, respectively, are obtained by summing over the polarization of the c, c quarks as, 

a-{t_L) = (y{tLCR) + a{tLCL), cr{t_R) = cr{t_BCL) + cr{t_RCR), 

a{tL) = cr{tLCR) + CT(iLCL), (j{tR) = cr{tRCL) + <7{tRCR). (14) 

Forming the half differences and sums of rates, 5a = ^{<y{tL) — o-{tR)), Sa — ■^{(j{iR) — a{tL)), (Jav — \{cr{tL) + 
(T(tij)), CTqu = \{<7{tR) -\- cr(ii))], such that, a{tL,R) = Uav ± 5a, a{tR^L) — o'av ± 5a, one can define the following two 
CP-odd combinations, 

^ ^ Tgy - agy ^ ate " ^fc j^poi ^ 5a -5a 

which will be designated as unpolarized and polarized integrated rate asymmetries. The above definition for the 
unpolarized asymmetry, A, is identical to the one studied in our previous work. M The asymmetries depend on 

the RPV coupling constants through the ratio of loop to tree amplitudes as, Im{ ,u"° ,]''''' ) ex: smip, where the 

IJfc' Ij'k' 

dependence on the CP violation angle parameter, f/;, reflects the particular prescription adopted in this study to 
include the CP-odd phase. The index k' refers to the d-squark generation in the tree amplitude and the indices i, k 
to the fermion-sfermion generations for the internal fermion-sfermion pairs, (g/), (r* ), in the loop amplitude. 

It is important not to confuse the above analysis with that of the top-antitop pair production, /"/"*" — > tt, where a 
CP-odd asymmetry observable for a single final state may be defined in terms of the difference of helicity configurations, 
f (titi) — a{tRtR). A non vanishing value for the corresponding difference of polarized rates can only arise via tree- 
loop interference terms involving the absorptive part of the top quark electric dipole moment, Im{d^''). Mm One 
should note here that the one-loop contribution of the RPV A' interactions to I'm{dj'^) vanishes. Two closely related 
processes, which are amenable to an analogous treatment, are the bb quark pair |p8| and r"^ t~ lepton pair production. 
Double spin correlation observables for the latter reaction, l~l^ — > t^t^ , have been examined in a recent work. |39] 
We note that the RPV A interactions can give a non vanishing contribution to Im{d^''). 

The results for the rate asymmetries are displayed in Fig.p. The numerical results for the unpolarized case (window 
{A) in Fig.g) update the results presented in ref. [Q since the present calculation includes the contributions from 
the Lorentz covariant tensorial (cr^ij coupling which were ignored in our previous work. |y] The asymmetry for the 
polarized case (window {B) in Fig.ra) involves the difference of the spin helicity asymmetry in the total production 
cross sections for the CP mirror conjugate top and antitop mirror reactions. While this CP-odd polarized asymmetry 
is not directly observable, it enters as an important intermediate quantity in evaluating the measurable kinematic 
distributions of the top decay products dependent on the top spin. We have assumed all the relevant RPV coupling 
constants to be equal and set the CP-odd phase at sini/; = 1. The rapid change in slope for the m = 200 GeV 
case are due to the threshold effect from the imaginary part in the superpartner one-loop contributions, which set at 
•\/s = 400 GeV . Aside from this large discontinuous contribution, one sees that both asymmetries comprise another 
contribution which is nearly independent of in and increases smoothly with the initial energy. Both asymmetries, A 
and A^°^ , take values of order a few 10~^, reaching O(10~^) at the highest incident energies. 

The statistical uncertainties on the asymmetry may be evaluated in terms of the signal cross sections and the 
integrated lunfinosity by considering the approximate definition, 5 A = l/[C{atc + <^ic)\^ ■ Using the same input value 
for the l umin osity C = 100. fb^^ at the three cm energies, ^/s = [0.192, 0.5, 1.] TeV, along with the cut signal rates 
in Table pi A| , we obtain statistical errors on the asymmetries of order 0(10^^). These values lie nearly two order of 
magnitudes above the value obtained for the signal. At this point, it is important to observe that in getting the above 
estimates for the rates we have been using somewhat conservative assignments for the RPV coupling constants. As 
already noted, the single top production cross sections could possibly be two order of magnitudes larger if we were 
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to use coupling constants values of order, A'j^jfe'^isfc — 10~^- Such values are compatible with the indirect bounds 
only for the extreme down squark mass in = 0(1 TeV) range. In the hypothetical case where the production rates 
would be enhanced by two order of magnitudes, the statistical errors on the asymmetries would correspondingly get 
reduced by a factor O(10~^), thereby reaching the same order of magnitude as the signal asymmetries. Nevertheless, 
as plotted in window (A) of Fig.O, the corresponding errors would still be somewhat larger than the signals. We 
should note here that the contribution to the one-loop amplitude from internal sfermion and fermion lines belonging 
to the third generation is controlled by the coupling constants quadratic product, A323A333, which is subject to weak 
constraints. Should the RPV coupling constants exhibit a hierarchical structure with respect to the quarks and leptons 
generations, one cannot exclude the possibility of a factor 10 enhancement from the ratio, /'7i(A323A333/A]^3A']^33). 
Such an order of magnitude gain on this ratio would raise the asymmetries up to O(10~^) bringing them well above the 
experimental uncertainties. Lastly, we observe that a more complete formula for the uncertainties on the asymmetries 
reads, {6A)^ = 2(^crfe)^[l — C+(l + C)^^]/(crfc + o'tc)^j where we used equal standard deviations for the CP conjugate 
reactions rates, Sate = Sate, and denoted the correlated error on these two rates as, C =< SatcSaic > /Sa^g. Clearly, 
an improvement on the statistical treatment of the tc + ic events sample, allowing for a positive non vanishing value 
of the error correlation associated with the identification of isolated single negatively and positively charged lepton 
events, should greatly help in reducing the experimental uncertainties caused by the small event rates. 
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FIG. 6. The CP-odd production rate asymmetries as a function of the center of mass energy, S2, for fixed values 
of the down squark mass, rh = [100, 200, 500, 1000] GeV. The left hand plot (A) gives the unpolarized asymmetry, 
A — {cFav — <Jav)/io'av + o^av)- The Upper bouuds for the absolute values of the statistical errors on the asymmetries, as 
evaluated with A'j^2fc-'^i3fc = 0.1, rh — 100 GeV and integrated luminosities £ — 100. fb~^ are shown as full circles. The right 
hand plot (B) gives the spin polarization dependent asymmetry, A^°^ = {5a — Sa)/{aav + aav)- 

The energy distribution for the negatively and positively charged leptons in the pair of CP-conjugate reactions may 
be defined as, 



< a' 



da+ 
>^< -fET >=<critL)fL + a{tR)fR>, <a >=< 
abji 



da 



-^ >=< a{tR)fL + a{tL)fR >, 



(16) 



where the correlations between the top spin lepton momentum are described by the factors, /l,_r = 2(1 "F cos0^), 
and the brackets stand for the angular integration. The occurrence of angular correlation factors of opposite signs in 
the t production case accounts for the kinematical fact that the antitop is oriented in space with a momentunr ~p. 
A CP-odd charge asymmetry observable with respect to the charged lepton energy distribution may be defined by 
considering the following normalized difference of distributions. 



^jpol ^ 



<a^ > - <a > 



{(Tav 



.)+ < (da ^ Sa)ifL - fR.) > 



< a+ > + < a- > [aav + crav)+ < [Sa + 6a){fL ~ Ir) >' 



(17) 



The numerical results for the charged lepton energy distributions and for the above defined charge asymmetry in 
the lepton energy distributions are displayed in Fig.^. (Note that the transverse energy distribution, in the plane 
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orthogonal with respect to the top momentum, may be simply obtained as, 



dT 



dr 



The distribution in the 



dEiT ~ dEi sin 6*^ ' 

plane orthogonal to the collision axis is less trivial to evaluate since this requires an additional integration over the 
lepton azimuthal angle.) The energy distributions for the unpolarized cross section essentially reproduce the results 
found in our above quoted event generator predictions, Fig.0. The energy distributions for the polarized asymmetry 
lie at values of order of magnitude, O(10~'^), always retaining the same positive sign as the lepton energy varies. For a 
fixed energy of the emitted lepton, the asymmetry increases with the initial energy, reaching values of order 0(10^^). 
In window {B) of Fig.[7| we have plotted the experimental uncertainties using the same inputs for the luminosities and 
the rates as in the discussion of the unpolarized asymmetries given above. To ease the comparison with experiment, 
we divide the charged leptons energy interval into three bins of width IQQGeV each, centered at the three lepton 
energies, Ei = (50, 150, 250) GeV . The statistical errors on the asymmetries in the energy distributions lie at the 
same level as those associated to the total asymmetries, so that similar conclusions should apply. Setting ourselves 
within the same optimistic scenario by using A'j^jfc-^isfc = 10~^ and C — 100/6""'^, we obtain expected errors of order 
O(10~^). These values are insufhcient for a comfortable identification of a signal asymmetry. However, we reiterate, 
as in the above discussion, that an enhancement of the signal asymmetries to an observable level of 0(10""'^), due to 
a hierarchical structure in the generation dependence of the A^ j., is a real possibility. 
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FIG. 7. Energy distribution for the charged lepton as a function of the laboratory frame lepton energy, for a set of center of 
mass energy, s^ = [192, 500, 1000] GeV. The parameters are set at. A' = 0.1, m = 100 GeV. The left hand plot (A) gives 
the differential lepton energy distribution, -j^. The right hand plot (B) gives the asymmetry in the energy distribution for 

leptons of opposite charges in the CP-conjugate final state channels, (tc) and (ci): AA''°'' = [^ ^~]/[^~ + im~^- "^^^ 

upper bounds for the absolute values of the statistical errors on the asymmetries, as evaluated with Xi2k^i3k = 0-1 S'^d with 
integrated luminosities, jC = 100. fb~^, are shown for three energy bins of width 100 GeV each, centered at the charged lepton 
energies, Ei = (50, 150, 250) GeV. The results for three values of the center of mass energy, S2 = [192, 500, 1000] GeV are 
displayed by full triangles, squares and circles. 

IV. CONCLUSIONS 

We have demonstrated that single top production through the RPV interactions could be observed at the future 
linear colliders or else be used to set bounds on the RPV coupling constants, yi2k^'i3k < O(I0~^), over a wide interval 
for the down squark mass, mj < 1. TeV. The b quark tagging would help greatly to overcome the background. Even 
with an imperfect b quark tagging, it is still possible to drastically reduce the background, from WW and 66, without 
much harming the signal. The analysis of top polarization observables via the semileptonic decay channel of the top 
allows to test for the presence of a CP violating complex phase, embedded in quadratic products of the RPV coupling 
constants. We have focused on the asymmetry in the energy distributions of the charged leptons in the CP-conjugate 
pair of final states, bl^vc and bl^Dc, obtaining asymmetries of order 10"'^ — 10~^ for the incident energies expected 
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at the future leptonic colliders. These values lie somewhat below the anticipated limits of observability. However, 
it may be possible to obtain enhanced values of order 10~^, should the RPV coupling constants A^ j. exhibit large 
hierarchies with respect to the quarks or leptons generations. Future promising extensions might include analogous 
reactions accessible with lepton-photon or photon-photon colliding beams, ^7 — *■ ic, 77 — > tc, where the expected 
production rates are substantially larger than those for the l^l^ colliders. 



APPENDIX A: 



Helicity amplitudes. 

The helicity spin basis Dirac spinors for a fermion or an antifermion, of mass m and four momentum, A:^ — 
{Ek — (fc^ + m^)2^ fc), and polar coordinates, k = {9, 0), can be written in the form of direct products of the Dirac 
spinor two-component space with the the two-component space of Pauli helicity basis spinors, (j)x{k), satisfying, 

r 



f? • k(j)\{k) — X(f>\{k). In the Dirac representation for the Dirac matrices, 70 
read, 



/3, 7 = f3d, 75 



-1 



, the spinors 



u{k, A) 
0-1 (fc) 



1 



fcA 

sin(6'/2)e-*"^ 
cos(6l/2) 



X 0A(fc), v{k,X) 
hiik) = 



-fcA 
1 

cos(6'/2) 
sin(0/2)e+*"^ 



X (f>-\ik), 



(A.l) 



where, €k ^ Ef^ + m, k = \k\/{Ei^ + m), and xai [-^ — ^1]; ^^^ the Pauli spinors in the basis with a fixed quantization 
axis identified with the spatial three-axis, Oz. The helicity basis spin eigenstates with a space parity reversed three 
momentum are defined as, 0a(-A;) = e-'('^+'')(-^'~^)/2(g-i^^<T„)^,^^^, ^ ,^A(fc)|[0^7r-e, 0^0+7^]- 

The 8 non vanishing helicity amplitudes for the process, l^{k' , /i') -I- l^{k, /i) —f uj(p, A) -t- uji{p' , A'), are listed in 
the formulas below: 



-+ 



= AT [-((1 + pp') (Xi + X2) + {p + p') (X3 + X4)] sin2(0/2), 
= 2T[{-l+pp'){X^+X2) + {p-p')X^ + {p-p')X^ 
+ 2p{p + p') (X^+Xe) + 2p{\+pp') (Xt + Xs)] sin(0), 
= -AT[~-{l+pp' ) (-Xi +X2) -{p + p') (X3 - X4)]cos2(0/2), 
= 2T [{-l+pp' ) {-X, + X2) + i-p + p' ) (X3 - X4) 
-2p{p + p') (X5 - Xe) - 2p{\+pp') (X7 - Xg] sin(0), 
= 2T [(-1 + pp') (Xi + X2) + {-p + p')X-i + i-p + p')Xi 
+ 2p{p + p') {X^+X^) - 2p{l+pp') (X7 - Xs)] sin(0), 
= -AT [(1 + pp') (Xi + X2) + {p + p') (X3 + X4)] cos2(0/2), 
= 2.F [{-l+pp' ) {-Xi + X2) + {p-p') (X3 - X4) 
- 2p (p + p') (X5 - Xe) + 2^(1+ pp') {Xj - Xs)] sin(0), 
Afg = Af(- + -+) = 4.F[-(1 +pp') (Xi - X2) - {p + p') {X3 - Xi)] sin2(0/2). 



Ml = M(+ - 
Afa = M{+ - 

A'h = M{- + 

Mi = M{- + 

M5 = M{+ - 



Me = M{+ - -+ 

Afr = M{- + — 



(A.2) 



The arguments refer to the fermions helicity in the following order, Mi{{hf,+ ,hf,- ,hf,hf). The remaining helicity 
amplitudes, omitted from the above list, are understood to vanish identically. We denote by 9 the top scattering 
angle, cos9 = k ■ p, by [Ep, E'^^ - ^. -u ^2 -, ^2 



(s ± m^j =F mj,)/2\/s, the top and charm quarks energies, and use the following 



abbreviated notations, p ~ -^-^ — 
notations, 



;;' _ P 



p = 



T = \\s{Ep 



mj 



){K 



nij 



)] 2 , along with the useful compact 



Xi = Ga+A + 7^, X2 = Ga-A + 7^, X3 = Ga+B + 7^, X4 = 

X5 = -Ga+C, Xq = -Ga-C, X-j = ica+P, Xg = ^GaTV, 

2 '2 2 '2 

where Ga^A, ■ ■ ■ are defined in eq. @, IZ in eq. (|^), and A,- ■ ■ ,V, in eq. (|^). 
One-loop RPV vector boson vertex functions. 

The one-loop vertex functions, as derived in W, are given by the formulas. 



Ga^B + Tl, 



A't' = 



^'iJ'k\*jk I 

(47r)2 I 



(A.3) 



,(2) 



aL{u)B['^+a{fL)mf^Ca + a{f'){2C2i + 2m'j{Ci2 - C21 + C23 - Cn )) 
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+ aifn) (4'^ - 2 C24 - mp^Co + m/{Co + SCn - 2Ci2 + 2C21 - 2C23) - m^Cu) 



A 



JJ' 



^' iJ'k\*jk 



-,JJ' 



mjmj, [2a(]') (-C23 + C22) + a(/fl) (-Cn + C12 - 2C23 + 2C22) ] 
± mj[a{fR){Cii - C12 + C21 - C23) 



(47r)2 



(47r)2 4 

a(/')(Cii + C21 - C12 - C2z)] + mj,[a{.fR){C22 - C23) + a(/')(C'23 - C22)] 



(A.4) 



The relevant configurations for the internal ferniion and sfermion propagating in the loop are: (^,) = {-*\ (/' )• 

The notations for the Passarino-Veltman two-point and three-point integrals, as specified in our work, [Q are defined 

according to the following conventions, S^ ~ BA{—p — p',Tnf,mf), B)^ = BAi—pTTnf^mi,)^ [A — 0,1] and 

Ca = CA{—p,—p',iTif,mf,,mf),CA = CA{—p,~p','mf,,mf,mr,). [A = 0,11,12,21,22,23] The integral functions 
with a tilde are associated, with the one-loop diagram for the sfermion current. 
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